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Abstract 



Let © n denote the symmetric group with n letters, and g(n) the max- 
imal order of an element of (S„. If the standard factorization of M into 
primes is M = q^q* 2 ...q™ k ,we define £(M) to be q™ 1 + q™ 2 + . . . + q™ k ; 
one century ago, E. Landau proved that g(n) = max(( M j<„ M and that, 
when n goes to infinity, log<;(n) ~ \J n log(n). 

There exists a basic algorithm to compute g(n) for 1 < n < JV; its 

running time is O ^7V 3//2 /\/log and the needed memory is O(N); it 
allows computing g(n) up to, say, one million. We describe an algorithm 
to calculate g(n) for n up to 10 15 . The main idea is to use the so-called 
i-superchampion numbers. Similar numbers, the superior highly composite 
numbers, were introduced by S. Ramanujan to study large values of the 
divisor function r(n) = ~^2 d , n 1. 

Key words: arithmetical function, symmetric group, maximal order, highly 
composite number. 

2000 Mathematics Subject Classification: 11Y70, 11N25. 

1 Introduction 

1.1 Known results about Landau's function 

For n > 1, let 6 n denote the symmetric group with n letters. The order of a 
permutation of 6 n is the least common multiple of the lengths of its cycles. Let 
us call g(n) the maximal order of an clement of & n . If the standard factorization 
of M into primes is M = q^q^ 2 ■ ■ ■ q^ k , we define £(M) to be 

£(M) = q^+q« 2 + ... + q^. (1.1) 

E. Landau proved in [H] that 

g(n) = max M (1.2) 



*Research partially supported by INRIA and by CNRS. 



1 



which implies 

£(g(n)) < n (1.3) 

and for all positive integers n, M 

i(M) < n =>■ M < g(n) <S=J> M > g(n) =>■ £(M) > n. (1.4) 
P. Erdos and P. Turan proved in [S] that 

M is the order of some element of & n £(M) < n. (1-5) 

E. Landau also proved in [9] that 



\ogg(n) ~ ynlogri, n — * oo. (1-6) 

This asymptotic estimate was improved by S. M. Shah and M. Szalay [3D]; 
in [12] , it is shown that 

logp(n) = ^Li _1 (n) + 0(y/nexp(-ay/logn)) (1.7) 

for some a > 0; Li -1 denotes the inverse function of the integral logarithm. 

The survey paper [TJj of W. Miller is a nice introduction to g(n); it contains 
elegant and simple proofs of (|1.2j) . (|1.5p and (|1.6p . 

J. -P. Massias proved in [TT] that for n > 1 

log 5 (n) < logg(131936 6) w i. 5313V^i^. (1.8) 

V ;_ ^1319366 1og(1319366) V v b v ; 

In [13] more accurate effective results are given, including 



logg(n) > \J n log n, n > 906 (1.9) 

and 

/ x / — ; / log log n — 0.975 \ 

Iogff(») < V^logn 1 + -2-f- , n>4. 1.10 

V 21ogn y 

Let P + (g(n)) denote the greatest prime factor of g(n). In [5J, J. Grantham 

proved 

P + (g(n)) < 1.328a/" log n, n > 5. (1.11) 

Some other functions similar to g(n) were studied in [7], [TD], [52], [3D] and [3T] , 



1.2 Computing Landau's function 

A table of Landau's function up to 300 is given at the end of [IS]. It has 
been computed with the algorithm described and used in [19j to compute g(n) 
up to 8000. By using similar algorithms, a table up to 32000 is given in [T5"] . 
and a table up to 500000 is mentioned in [8]. The algorithm given in [19] will 
be referred in this paper as the basic algorithm. We shall recall it in Section 
[21 It can be used to compute g(n) for n up to, say, one million, eventually a 
little more. It cannot compute g(n) without calculating simultaneously g(n') 
for 1 < n' < n. 
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If we look at a table of g(n) for 31000 < n < 31999 (such a table can be 
easily built by using the Maple procedure given in Section [2|), we observe three 
parts among the prime divisors of g(n). More precisely, let us set 

g(n) =Y[p a », gM{n) = J[ p<*», g™(n) = J[ p"*, g^(n) = J[ p a »; 

p p<17 19<p<509 p>509 

the middle part g^ 2 \n) is constant (and equal to ni9<p<509^) f° r a ^ n between 
31000 and 31999, while the first part g^ 1 ' (n) takes only 18 values, and the third 
part (n) takes 92 values. 

So, if n' is in the neighbourhood of n, g{n') / g(n) is a fraction which is the 
product of a prefix (made of small primes) and a suffix (made of large primes) . 

The aim of this article is to make precise this remark to get an algorithm 
able to compute g(n) for some fixed n up to 10 15 . 

1.3 The new algorithm 

Let r(n) = | n 1 be the divisor function. To study highly composite numbers 
(that is the n's such that m < n implies r(m) < i~(n)), S. Ramanujan (cf. [2"4l 
I25[ I20|) has introduced the superior highly composite numbers which maximize 
r(n)/n E for some e > 0. This definition can be extended to function t: N is said 
to be €-supcrchampion if it minimizes £(N) — plog(A r ) for some p > 0. These 
numbers will be discussed in Section [5] they are easy to compute and have the 
property that, if n ~ £(N), then g(n) = N. 

If N minimizes £(N) — p\og(N), we call benefit of an integer M the non- 
negative quantity ben (M) = £(M) — £(N) — p\og(M/N). If n is not too far 
from i(N), a relatively small bound can be obtained for ben c/(n), and this allows 
computing it. This notion of benefit will be discussed in Section [5] 

To compute g{n), the main steps of our algorithm are 

1 . Determine the two consecutive €-superchampion numbers N and N' such 
that £(N) < n < £(N') and their common parameter p (cf. Section [5|). 

2. For a guessed value B' , determine a set V(B') of plain prefixes whose 
benefit is smaller than B' (cf. Section [7Tl and Section 17721 . 

3. Use the set T>(B') to compute an upper bound B such that ben g(n) < 
beng(n) + n - £(g(n)) < B (cf. Section [773]): note that, from ([Q]) . 
£(g( n )) < n holds. 

4. Determine T>(B) : a set containing the plain prefix of g{n). If B < B' , to 
get T>(B), we just have to remove from V{B') the elements whose benefit 
is bigger than B. If B > B', we start again the algorithm described in 
Section [7721 to get T>(B') with a new value of B' greater than B. 

5. Compute a set containing the normalized prefix of g(n) (cf. Sections 17. 7\ 
17781 and 1779)1 . 

6. Determine the suffix of g(n) by using the function G(pk, m) introduced in 
Section 11.41 and discussed in Sections [8] and [9] 
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In the sequel of our article, " step " will refer to one of the above six steps, and 
" the algorithm " will refer to the algorithm sketched in Section 11.31 

On the web site of the second author, there is a Maple code of this algorithm 
where each instruction is explained according with the notation of this article. 

If we want to calculate g(n) for consecutive values n = n\, n = ni + l, . . . , n = 
n%, most of the operations of the algorithm are similar and can be put in com- 
mon; however, due to some technical questions, it is more difficult to treat this 
problem, and here, we shall restrict ourselves to the computation of gin) for one 
value of n. 

To compute the first 5000 highly composite numbers, G. Robin (cf. [ST] ) 
already used a notion of benefit similar to that introduced in this article. 

1.4 The function G(pk, m) 

In step 6, the computation of the suffix of g(n) leads to the function G(pk,m,), 
defined by 

Definition 1. Let pk be the k-th prime, for some k > 3 and m an integer 
satisfying < m < Pk+i — 3. We define 

t~<l \ QlQ2 ■ Qs / -. -1 rj \ 

G(pk,m) — max (1-12) 

gi<?2 • • -q s 

where the maximum is taken over the primes Q\, Q2, • ■ • , Q s , qi, 12, ■ • • , Qs (s > 
) satisfying 

3 < q s < q s -i < ■ ■ ■ < qi < Pk < Pk+i < Qi < Qi < ■ ■ ■ < Qs (1-13) 

and 

s 

1=1 

This function G(pk,m) is interesting in itself. It satisfies 

e{G( Pk ,m)) <m. (1.15) 

We study it in Section [SJ where a combinatorial algorithm is given to compute 
its value when m is not too large. For m large, a better algorithm is given in 
Section El 

Let us denote by ni(n) < H2{n) < . . . the increasing sequence of the primes 
which do not divide g(n), and by P(n) the largest prime factor of g(n). It is 
shown in |17j that limn.—^ P(n)/fj,i (n) = 1. We may guess from Proposition 
[inithat /xi(n) can be much smaller than P(n) while U2(n) is closer to P(n). It 
seems difficult to prove any result in this direction. 



1.5 The running time 

Though we have the feeling that the algorithm presented in this paper (and 
implemented in Maple) yields the value of g(n) for all n's up to 10 15 (and 
eventually for greater n's) in a reasonable time, it is not proved to do so. 

Indeed, we do not know how to get an effective upper bound for the benefit 
of g(n) (see sections El 17.31 and lll.lj) and in the second and third steps, what 
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we do is just, for a given n, to provide such an upper bound B = B(n) by an 
experimental way. 

In the fourth step, the algorithm determines a set ~D{B) of plain prefixes 
(cf. sections 17.21 and l7.3p . It turns out that the number v{n) of these prefixes 
is rather small and experimentally satisfies v(n) = O(n - 3 ) (cf. (|7.11[) ); but we 
do not know how to prove such a result, and it might exist some values of n for 
which v{n) is much larger. 

Let us now analyze each of the six steps described in Scction [1.3l 

The first step determines the greatest supcrchampion number ./V such that 
< n. Let S(x) = ^2 p<x p be the sum of the primes up to x. The main 
part of this step is to compute S(x) for x close to y/nlogn. In our Maple 
program, by Eratosthenes' sieve, we have precomputed a function close to S(x), 
the details are given in Section [5l However, a faster way exists to evaluate 
S(x). By extending Meissel's technique to compute n(x) = J2 P < X 1> ( CI - 02)> 
M. Deleglise is able to compute J2 P < X f(p) where / is a multiplicative function. 
E. Bach (cf. [TJ [5]) has considered a wider class of functions for which this 
method also works. By his algorithm, M. Deleglise has computed S^IO 18 ), and 
S(x) costs 0(x 2 / 3 / log 2 x). We hope to implement soon this new evaluation of 
S(x) in our first step. 

The second and the fourth steps compute respectively V{B') and T>{B). If B' 
is "well" chosen, we may hope that Card(2?(B')) is not much larger than v(n) = 
Card(I?(_B)). The running time of the computation of T>(B') as explained in 
Section 1731 could be larger than v(n). For n m 10 20 , most of the time of the 
computation of g(n) is spent in the second and fourth steps. But any precise 
estimation of these steps seems unaccessible. 

The running time of the third step is 0(Card(P(_B'))), and we may hope 
that it is 0(v(n)). 

In practice, the fifth step (finding the possible normalized prefixes) is fast. 
For every plain prefix tt, Inequations (|7.36|) have at most one solution, and the 
cost of this step is 0(v(n)). 

The sixth and last step also is fast. Under the strong assumption that 5\{p) 
is polynomial in logp (see (|9.8p ). for any m, the computation of G(jp, to) (where 
p is a prime satisfying p m \J n log n) is polynomial in logn, and the number of 
normalized prefixes surviving the fight (cf. Section 17. 9|) seems to be bounded 
(we have no examples of more than three of them), so that (see Section [778]) this 
step might be polynomial in logn. 

1.6 Plan of the paper 

In Section [3J some mathematical lemmas are given. The various steps of the 
algorithm presented in Section 11.31 are explained in sections HE] Section fTUl 
presents some results while Section [TT] asks five open problems. 

1.7 Notation 

We denote by V = {2, 3, 5,7,.. .} the set of primes, by p G V a generic prime, 
by pi the z-th prime and by v p (N) the p-adic valuation of N, that is the greatest 
integer a such that p a divides N. Qi and qi also denote primes, except in Lemma 
[T]which is stated in a more general form, but which is used with Qi and qi primes. 
The integral part of a real number t is denoted by \t\ . The additive function £ 
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can be easily extended to a rational number by setting £(A/B) = 1(A) — 1(B) 
(with A and B coprime). 



2 The basic algorithm 
2.1 The first version 

For j > 0, let us denote by Sj the set of numbers having only pi,P2, ■ ■ ■ ,Pj as 
prime divisors 

Sj -.={M;p\M => p<p 3 }. (2.1) 

We have S = {1}, Si = {1,2,4,8,16,...}. The algorithm described in [15] 
computes the functions 

qAn) = max M (2.2) 

MeSj, £(M)<n 

which obviously satisfy the induction relation 

9j( n ) = max [g j -i(n),p j g j -- L (n - pj), . . . ,p$gj_i(n - p*)] (2.3) 

where k is the largest integer such that p^ < n, and go(n) = 1 for all n > 0. 
Using the upper bound (jl.llj) . we write the following Maple procedure: 

Algorithm 1 The basic algorithm: this Maple procedure computes g(n) for 
< n < N and stores the results in table g. 

gden:= proc(N) local n, g,pmax,p,k,a 

for n from to N do 
g[n] := 1 

endo; 

pmax := //oor(1.328*etiaZ(sgrf(iV*logiV))); 
p:=2; 

while p < pmax do 

for 7i from N to p by — 1 do 
for k from 1 while p k < n do 

a := p k -k g[n — p k ]; 
if g[n] < a then 

g[n) := a 
end if 
endo 
endo; 

p: =nextprime (p) 
end while; 
end: 



The running time of this procedure is 13 hours for N = 10 6 on a 3 Ghz 
Pentium 4 with a storage of 337 Mo. To compute g(n), 1 < n < N, the 
theoretical running time is O (iV 3 / 2 / \J\og N) and the needed memory is 0(N) 
integers of size exp(0(y / iVToglV)). 
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2.2 The merging and pruning algorithm 

The above algorithm takes a very long time to compute gj (n) when j is small. It 
is better to represent (gj(n)) n >i by a list Lj = [[Mi, li], . . . , [M i: k], . . .] (where 
U = £(Ali)) ordered so that Mi + i > Mi and U+\ > k- If U < n < k+i, then 
9j (n) = Mi. So, Lq = [[1, 0]] and L 1 = [[1, 0], [2, 2], [4, 4], [8, 8], . . .]. 

To calculate Lj+i from Lj we construct the list of all elements [MiPj +1 , li + 
£(Pj +1 )] for all elements [M % ,k] € Lj and a > such that k + £(Pj +t ) < N. 
We sort this new list with respect to the first term of the elements (merge sort 
is here specially recommended) to get a list A = [[Ki, Ai], [K%, Aa], • • •] with 
K\ < K<i < . . . Now, to take (|2.3p into account, we have to prune the list A: if 
K r < K s and A r > A s , we take off the element [K r , X r ] from the list A. The list 
Lj + i will be the pruned list of A. 



3 Two lemmas 

Lemma 1. Let s be a non-negative integer, and t\, qi, q2, ■ ■ ■ , q s , Qi, Q2, ■ ■ ■ , Qs 
be real numbers satisfying 

< h < q s < q s -! < . . . < qi < Qi < Q 2 < ■ ■ ■ < Q s - (3.1) 

s 

If we set S = Qi — qi, then the following inequality holds: 

i=l 

, Q1Q2...QS . (s\ 

1. < cxp — . 

q\q2---q s \t\) 

Moreover, if s > 1 and S < Q±, we have 

Q1Q2 ■ ■ ■ Q s , Qs Qs-i Qi 
■ ~~ ~~ *c ~~ — < c. . . . ^ 



qiq 2 ...q s Qs - S Q s -i - S Ql - S 

with the first inequality in 2. strict when s > 2. 

Proof. Lemma[T]is a slight improvement of Lemma 3 of |18j where, in 2., only the 
upper bound Q\ j (Qi — S) was given. Point 1. is easy by applying 1 + u < exp u 
to u = Qi/q% — 1. Let us prove 2. by induction. For s = 1, 2. is an equality. 
Let us assume that s > 2. Setting S' = J2t=2 Qi ~ li = & — (Q 1 — we nave 
5' < S < Qi < Q s and by induction hypothesis, we get 

Q1Q2 ■ ■ - Qs _ Qi Q2 ■ ■ ■ Qs . Qi Qs /„ „n 

qi 92 ■ ■ • q s qi qi---q s qi Qs - & 

We shall use the following principle: 

Principle 1. If x andy add to a constant, the product xy decreases when \y—x\ 
increases. 

We have Q s — S' < Q s — (Q s — qs) = qs < qi, and using Principle 1, we get 
by increasing qi to Q\ and decreasing Q s — S' to Q s — S 

qi(Q s -S')>Q 1 (Q s -S) 

which, from (|3.2|) . proves 2.. □ 
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y — y 

Lemma 2. Let x > 4 and y = y(x) be defined by = The function 

logy logx 

y is an increasing function satisfying y(x) > 2 and 

flf log 2 (4 + log2)log2 / 1 \\ 

1. y[x) = \ — 1 — — o h O — 5 — , x — > oo 

yw V 2 V 21oga; 81og 2 x Vlog 3 x/7 

U. y(.T) < /or x > A. 
3- y{x) <^forx> 80. 
Proof. 1. and 3. are proved in [T^], p. 227. Since i (f 2 — t)/ log i is increasing 

X — \/ X X 

for i > 1, in order to show 2., one should prove > which holds for 

5 log x log x 

x > 4. □ 



4 The superchampion numbers 

Definition 2. An integer N is said I- superchampion (or more simply super- 
champion) if there exists p > such that, for all M > 1 

e{M)-plogM>£(N)-plogN. (4.1) 

When this is the case, we say that N is a I- superchampion associated to p. 

Geometrically, if we represent logM in abscissa and £(M) in ordinate, the 
straight line of slope p going through the point (log M,£(M)) has an interscp 
equal to £(M) — plog(M) and so, the superchampion numbers are the vertices 
of the convex envelop of all these points (see Fig. [1]). 

Similar numbers, the so-called superior highly composite numbers were first 
introduced by S. Ramanujan (cf. [24]). The ^-superchampion numbers were 
already used in [13 H3 HH [12 HH [^] . The first ones are (with, in the third 
column, the corresponding values of p) shown in Fig. [2] 

Lemma 3. If N is an £- superchampion, the following property holds: 

N = g(£(N)). (4.2) 

Proof. Indeed, let N be any positive number and n = £(N); it follows from (|1.4p 
that N < g(n) — g(£(N)). If moreover N is a ^-superchampion, then, for all M 
such that £{M) < n = £(N), from P~T|) . we have plogM < plogN + £(M) - 
t(N) < p log TV which implies M < N, and thus, from (fl~2|) . (j4~2]) holds. □ 

Definition 3. 

1. For each prime p G V , let us define the sets 

£'={jP-\ S» = ii-lP pl+1 ~ pl \ S=£'U8" 



(4.3) 
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2 4 6 8 10 

log(N) 



Figure 1: The points (log(iV), t(N)), with C(N) < 50, for 1 < N < 60060. 



2. And we define 



£'=\J £' p , £" = (J £p and £=£'(J £" . (4.4) 
pev pev 



Remark: Note that all the elements of £ p are distinct at the exception, for 
2 2 2 - 2 

p = 2, of = — and that, for p ^ q, £ v D £„ = holds. 

log 2 log 2 

Lemma 4. Let p a real number. 

1. If p 6 £ p) p ^ I3§2 j there exist exactly 2 superchampion numbers associated 
to p. Let be N p the smaller one and the bigger one. Then iV+ = pN p 
and 



1 t P ^ ^ P 2 -P 



N p = p ap with a p = < 

P/ logp <p 



N+ = p a p with a+ = 

p/ logp<p 



logp 



logp 



V < P < 



logp 



logp 



(4.5) 



., P , ^ P -P 

if 7—<P< 



logp 



logp 



z / — ^ P < 



\ogp 



logp 
(4.6) 
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2 2 2 - 2 

2. If p = = € £, there exists superchampion numbers associated 

log 2 log 2 

to p: N p defined by (|4.5j) is egua? to 3, iV+ defined by (|4.6p is egitaZ to 12 
and i/*e i/iirci one is 6. 

5. //p £, i/iere exists a unique superchampion number N p = iV+ associated 
to p. Its value is given by both formulas ()4.5|) and ()4.6p . Lei // and p" be 
the two consecutive elements of £ such that p' < p < p" . Then we have 
N P = N P „=N+. 

4- Let us consider the sequence p^' defined by p^ = — oo, p^ = 3/log3 ; 
pW = 2/ log 2, p^ = (2 2 -2 1 )/log2 = pW, p<4> = 5/log5 and such that 
{p {i \i > 1} = £ and p« > pO" 1 ) /or i > 4. J/ A^°> = 1, 7V« = 3 ; 

= 6, A^ 3 ) = 12, A^ 4 -* = 60, etc... is the increasing sequence of all 
superchampion numbers, it satisfies: 

(i) For i > 0, N® divides iV(' i+1 ) and the quotient N^+V /N® is a 
prime number. The number of prime factors of , counting them 
with multiplicity, is equal to i. 

(ii) For i ^ 2, we have = N + {i) = N m+i) where N + (z) and iV_{<+i) 

are defined respectively in (|4. 5[) and (|4.6[) . 
(Hi) For all i > 0, is associated to p if and only if p^ < p < p' t+1 ) . 

(iv) Ifi^l (i.e., iVW ^ 3;, then v p (N^) is a non-increasing function 
of the prime p. 

Proof. We are looking for an A~ = nP Qp which minimizes F(N) = £(N) — 
p\ogN. 

An arithmetic function h is said additive if h{MxM?) = h{M\) + h(M2) 
when M\ and M2 are coprime. The functions log and £ are additive. Thus 
F is additive, and to minimize F(N) = ^2 p | N F(p v *( N >) we have to minimize 
F(p a ) on a for each p <G V . We have F(l) = and for p prime and i > 1, 
F{p l ) = p l — pilogp. The difference 

p — plogp if i = 

F{p*+ l ) - F{tf) = { (4.7) 
p l {p — 1) — plogp if i > 

is a non decreasing function of i that tends to +00 with i. Thus if F(p) = 
F(p) — F(0) = p—p\ogp > 0, the smallest value of F(p a ) is obtained for a = 0. 
If F(p) < let i be the largest positive integer such that F(jp l ) — F(p 1 ^ 1 ) < 0. 
Then the smallest value of F(p a ) is obtained on the set {j < i \ F(p>) = F(p 1 )} 
and the number of choices for a p is the cardinal of this set. 

This proves that we have more than one choice for the exponent a p if and 
only if there exists i > such that F(p i ) = F(p i+1 ). Due to (|4.7[) this is the 
case if and only if p G £ p . Moreover, the sets £ p being disjoint, there exists at 
most one p for which there are more than one choice for a p . 

If p > 3 we have p < (p 2 — p) < (p 3 — p 2 ) < • • • and there is at most one i 
such that F(p l+1 ) — F(p l ) = 0, so there are at most two choices for a p . 

For p = 2 we have 2 = 2 2 - 2 < 2 3 - 2 2 < • • • and for p = 2/ log 2 we have 
F(l) = F(2) = F(2 2 ), so we can choose for an every one of the three values 
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60 


12 


5/ log 5 
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9 
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19 
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9 
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30 
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9 
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60060 


43 


13/ log 13 
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ps 5.46 



Figure 2: The first ^-superchampion numbers. 



0, 1, 2. With this value of p we have F(3) = 3- (2/ log 2) log 3 < and F{p) > 
for p > 5. Thus there are 3 superchampion numbers associated to p = 2/ log 2 
which are 3,6, 12. This proves 1., 2., 3. and 4.; for more details, see [TH] , 

□ 

Lemma 5. Let p satisfy p > 5/log5 ~ 3.11. There exists a unique decreasing 
sequence (xj) = (xj(p)) such that x\ > exp(l) and, for all j > 2, Xj satisfies 
xj > I and 

j 7 — 1 

arl- - arl- xi 

-1 — = i = P- (4-8) 

logXj logxi 

We /lave a?so 

xi > 5 and x 2 > 2. (4.9) 

Proof. The uniqueness of xi results from p > exp(l) and the fact that t <— > 
t/logt is an increasing bijection of [exp(l), +oo[. The uniqueness of Xj for 
j > 2 comes from the fact that t i— > — t J_1 )/logi = fc 7-1 ^ — l)/logf is an 
increasing bijection of ]1, +oo[. The inequality Xj > Xj+± for j > 2 comes from 
the increase of j i— ► (t J — f J_1 )/logi for each t > 1. 

Let us prove that xi > X2. The definition Q4.8[) of X2 implies 



x 2 2 2 2 - 2 

= p> r -— = - — — ps 2.89 



log x 2 log 2 log 2 

With the increase of t 1— > (t 2 — t)/logi this proves x 2 > 2. Thus x 2 — x 2 > x%, 
and therefore 

x 2 x\ — x 2 Xi 



log x 2 log x 2 log xi 

which, with the increase of £ > t/logt on [exp(l),+oo[ yields x 2 > xi and the 
decrease of (x„). Finally x\j logxi = p > 5/ log 5 gives x\ > 5. □ 

Proposition 1. Let p be a real number satisfying p > 5/log5 ; -/V p t/ie smallest 
superchampion number associated to p and iV+ t/ie largest superchampion num- 
ber associated to p (cf. Lemma^. Then, with Xj as introduced in Lemma\5^ 
we have 

n p = u n v* and n p=ii n p 3 - ( 4 - iq ) 

j>l Xj + 1 <p<Xj j>l X j + 1 <p<Xj 



11 



i 


T[i].q 




T[i].p 


T[i\.£ 


1 


2 


2 


3 


7 


2 


3 


2 


13 


49 


3 


2 


3 


13 


53 


4 


2 


4 


43 


301 


5 


5 


2 


47 


368 


6 


3 


3 


67 


626 


7 


7 


2 


97 


1160 


8 


2 


5 


107 


1487 


9 


11 


2 


251 


6307 


10 


2 


6 


251 


6339 


11 


3 


4 


271 


7453 



Figure 3: The first elements of table T associated to E 2 . 



Proof. Due to (|4.5|) . a p = 1 holds if and only we have 



P <P<^ (4.11) 



logp logp 
and by the definition (|4.8j) of xi and X2, this is equivalent to 



P ^ xi x\ - x 2 p 2 - p 
< and — < 



logp logxi logx 2 logp 

By the increase of 1 1— > t/ logt on [exp(l), +oo[ and 1 1— > (t 2 — t)/ \ogt on [1, +oo[, 
this proves that for p > exp(l), a p = 1 holds if and only if X2 < p < aci. It 
remains to prove that, when p = 2, this equivalence is still true. In this case, 
2/ log 2 = (4 — 2)/ log 2, and (14. lip is never satisfied. By (O 

we have x 2 > 2, 

and x 2 < 2 < x\ is false. Thus, for every prime p, we have a v = 1 if and only if 
£2 < p < X\. 

p % — p 1 p l+l — p l 

For i > 2, a p = i if and only if < p < , and, by the 

logp logp 

definition (|4.8[) of and Xi+\ this is equivalent to 

p p < Ej and < p p 



logp log a;, logccj+i logp 

or Xi+i < p < Xi. This proves the first equality (|4.10p . The second one can be 
proved by the same way. □ 

5 First step of the computation of g(n): getting 
p, N, N f . 

5.1 Fixing our notation 

When p = 5/log5 we have N p = 12 and £(N p ) = 7 (sec Fig. 2). 
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Definition 4. From now on, n > 7 i/ra/Z &e a fixed integer, and our purpose is 
to compute g(n). We will denote by p the unique real number p £ £ such that 
p > 5/ log 5 and 

£(2\g < n < £(iV+). (5.1) 
We will also fix the following notation. 

1. N = N p , N' = N+ and N = Y\_P a " is the standard factorization of N . 

p 

2. We define x\ = x\{p) > 5 and x 2 = x 2 (p) > 2 by (|4.8j) . 

3. Let pk be the largest prime factor of N = N p . It follows from (|4.10p that 

Pk < xi < Pk+i (5.2) 
and, actually, x\ = Pk+i unless p € £" (in this case p^ < x\ < pk+i). 
4- Let us define B\ by 

Si =min^_2a- 2 ,y-V^r) >0. (5.3) 

We have 

2 < x 2 < \fx~[ < p < x\. (5.4) 

Let us prove (|5 .4[) . Inequalities (|4.9p give 2 < x 2 - With Lemma Point 2., it 
yields x 2 < y/xi. Since for all t > 1, \ftj logt > e/2 > 1 we have \fx x j logxi > 1 
and thus p = xi/\ogx\ > \fx x . 



5.2 The superchampion algorithm 

Given n, as already said, the first step in our computation of g(n) is to calculate 
p, N, N', x\, x 2 ,Pk, B\ as introduced in Definition fij 

We begin by precomputing in increasing order the first elements of £" and 
stop when we get the first r e £" such that £(N+) > 10 15 . We get a set E 2 
with 1360 elements, 

E - f 22 ~ 2 32 ~ 3 23 ~ 2 ' 
2 _ \ log2 ' log3 ' log2 

We construct a table T, indexed from 1 to card(i^) = 1360. Let r = 
(q3 +1 - q*)/ logo the i th element of E 2 . Then T\i] = [q,j,p, I] where I = t{N+) 
and p is the largest prime p such that pj logp < r. The supcrchampions following 
are obtained by multiplying it successively by the primes following p. Figure 
Ogives the first values of T[i] . (In the Maple program the T[i] 's are the elements 
of the table listesuperchE2). 

The superchampions that are not of the form N+ for an r £ E 2 can easily be 
obtained from this table. For instance, the successive values of £(N) between 368 
and 626 are 368 + 53 = 421, 421 + 59 = 480, 480 + 61 = 541 and 541 + 67 = 608. 

Two elements of £ can be close. For instance, the smallest difference between 
two consecutive elements of £ less than 8 • 10 9 is 

43083996283 144589 2 - 144589 
log 43083996283 log 144589 

= 1759505912.7146899772- 1759505912.7146800938 = 0.0000098834 



13 



and thus, working with 20 decimal digits is enough to distinguish the elements 
of £. For any n up to 10 15 , Algorithm 2 below determines the superchampion 
N = N p as defined in DefintionJU 



Algorithm 2 : computes N = N p for a given n < 10 15 . 
Construct table T. 

i : = the largest index such that T[i].i < n. 
£' := T[i + 1].£, q' = T[i + l].q, f = T[i + 

{r' = (q'i'-q'V'-^/log q' is the smallest element in Ei such that £(N r r) > n} 
t :=£' -qftf-Vtf-iy, 

{This is the value £(N) of the superchampion N preceding A+} 
if t < n then 

p := r' 
else 

tiq:= T[i].£ + nextprime(T[i].p); 
while no < n do 

p := nextprime(p); hq := uq + p 
end while 
P :=p/logp 
end if 



6 Benefits 

6.1 Definition and properties 

Definition 5. Let p € £ and N = N p (as defined in Definition^. If M is a 
positive integer, from (|4.ip . we have £(M) — plogM > £(N) — plogN. We call 
benefit of M the non-negative quantity 

M 

ben (M) = £(M) - £{N) - plog — • (6.1) 
Let M = Yi p P^ p be the standard factorization of M. We define 

ben p (M) = £{p^) - £(p a ") - p{[3 v - a p ) logp > 0, (6.2) 

which implies 

ben (M ) = bcn P ( M )- ( 6 - 3 ) 

p 

Geometrically, if we represent logM in abscissa and £{M) in ordinate, the 
straight line of slope p going through the point (log M,£{M)) cuts the y axis at 
the ordinate yu = £(M) — p log(M) and so, the benefit is the difference yu — yn 

£(N') - £(N) 

(see Fig. H. Note that p = — y - — '- with N = N and N' = N+. 

v ^ log N' — log N 9 9 

Lemma 6. Let p g V . a = a p = v p (N) and 7 a non-negative integer. Then, 

1. ben(p 7 A) = £{p~ l+a ) — £{p a ) — pjlogp is non- decreasing for 7 > and 
tends to infinity with 7. 
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log(N)log(M) 



Figure 4: A = (log N,i(N)) and B = (log Af, £(M)). 



2. ben {N /p 1 ) = pj logp + £(p a 7 ) — ((p a ) is non- decreasing for < 7 < a. 
Proof. 1 . If 7 + a > 1 , we have 

ben (p~ f+1 N) — ben (p J N) = logp ( p 7 — : p 

V 1o sp 

which is non-negative from (|4.5p and tends to infinity with 7. 

If a = 7 = 0, we have ben (pN) — ben (iV) = log p(p/ logp — p) which is 
also non-negative from (|4.5|) . 

2. If a > 2 and < 7 < a — 2, we have 

/ N \ (N\ ( 1 p^-p*- 1 



vP 7+1 / VP 7 / V p 7 logp 

which is non- negative from (|4.5D . 
If a > 1 and 7 = a — 1, 

bcn ( 37+T ) ~ bcn ( 37 ) = lo S p 



vP 7+1 / VP 7 / V lo SP 

yields the same conclusion. 



□ 



Lemma 7. Let {//V 6e an irreducible fraction such that V divides N (as fixed 
in Definition^ and U = U-JJ 2 , V = V X V 2 with {U X ,U 2 ) = (Vi,V 2 ) = 1. Then 
we have 



1. 
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h ™{— J =bc H"^rJ +bcn lirJ- (6 - 5) 

Proof. Observing that a prime p divides at most one of the four numbers 
Ui,U 2 , Vi,V 2 we get ([6~4|) . By the additivity of the logarithm, ([63]) follows. □ 



The following proposition will be useful in the sequel. 

Proposition 2. Let M be a positive integer such that £{M) < n (thus, from 
CE3J, M < g(n) holds). Then, 

ben g(n) < ben M + £{g{n)) - £{M) 

and 

ben g(n) < ben g(n) +n- £{g{n)) < bcniU + n - £{M). (6.6) 
Proof. From (|6.1|) . we have 

ben 3 (n) - ben M = £(g(n)) - £{M) - plog < l{g{n)) - £(M) 

which implies the first inequality while the second one follows from (|1.3|) . □ 
We shall use Proposition [5] to determine an upper bound B such that 

ben g(n) < ben g(n) + n - £(g(n)) < B. (6.7) 
It has been proved in [13] that B < x\ and 

B = 0^=0( P ), (6.8) 

and, by the method of [23], it is possible to show that B = o(p). The largest 
quotient (ben g(n) +n — £(g(n)))/p that we have found up to n = 10 12 is 1.60153 
for n = 45055780. 

6.2 The benefit of large primes 

Proposition 3. Let N,B\,x\ and x 2 as in Definition^ If M is an integer 
satisfying ben (M) = £(M) - £(N) - p\og{M/N) < B x , we have 

1. if \fx\ < p then v p (M) < 1 

2- if x 2 < p < ^/xi then v p {M) < 2. 

Proof. 

1. Let us assume that the prime p satisfies p > ^/x 1 and divides M with 

exponent k > 2. With (|5.4|) . we have p > x 2 and, from (|4.10|) . the 
exponent a p of p in TV = N p is or 1. If a p = 1, from (|6.2p and (|4.5|) we 
have 



ben p (M) = p k — p — p(k — 1) logp = logp^^ 



p — p 

p 



i=2 

„2 



logp 



> l0 &p( P \ ogp P ~ P) =P 2 P l °SP ( 6 - 9 ) 
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while, if a p = 0, 

ben p M = p fc-pfclogp = logp^-p + g( g^ -p)) 

> i°gp f n — - - p) =p 2 -p- p logp. 

V logp / 

So, in both cases, (fO]) and (fo~2"l) yield ben M > ben p M > f(p) with 
f(t) = t 2 -t- plogt. We have /'(*) = 2t - 1 - p/t, /" (t) > and, as 
x 2 > 2 holds, (|4.8p implies 



f(x 2 ) = 2x 2 - 1 - ^— ^ > x 2 ( 2 - - J— ) 1 > 2 f 2 ) 
logx 2 V loga^/ V lo s 2 / 

and /(i) is increasing for t > x 2 . Thus, since p > y/xi, 



ben M > f(p) > /( y^i ) =x\- ^/x{ - - Xl log ^fx{ = — !- - sfxi > 5i 

logxi 2 

in contradiction with our hypothesis, and 1. is proved. 

2. Let p satisfy 2 < x 2 < p < \fx~i so that, from (|4.10p . a p = v p (N) = 1; let 
us assume that k = v p (M) > 3; one would have as in (|6.9p 

ben M > \ogp > ( — p I > p 3 — p 2 — p logp. 

^ V iogP / 

The function /(t) = t 3 — t 2 — plogt is easily shown to be increasing for 
t > x 2 . From (|4.8[) , /(x 2 ) = x\ — x\ — — x 2 ) and thus 

benM > x\ — .Tj — (x 2 — £2) == x 2 (x 2 — 2^2 + 1) > x 2 — 2x 2 . 



From (|5.3p , it follows that ben M > B\ holds, in contradiction with our 
hypothesis, and 2. is proved. 

□ 

7 Prefixes 

7.1 Plain prefixes and suffixes 
Definition 6. Let j be a positive integer. 

1. For every positive integer M let us define the fraction 

7T«(M) = [] p v P {M)-v p {N) = -Q p v P (M)-a p (? ^ 

P<Pi P<Pj 

and call w^'(M) the j -prefix of M. 

2. We note Tj, and call it the set of j-prefixes, the set of fractions 

T 3 = \ S = HP**! «p>-«p|- (7-2) 
I p<Pj ) 
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3. For B' > 0, we define 

T,j{B') = {5eTj; ben(NS) < B'} . (7.3) 

Definition 7. Le M be a positive integer. Let us define 

tt(M)= Yl p»^ M )- a » = 7 r^' l) (M) (7.4) 
P<V^T 

where Pj 1 is the largest prime less than yjx~i, and £(M) = M/(Nn(M)). Thus 
we have 

M = Nn(M) f(M). (7.5) 
7r(M) wzH fee called the plain prefix of M , and £(M) the suffix of M . 

Let us show that, for each j such that < >Jxi, we have 

ben (iV7r (1) (A/)) < . . . < ben (Ntt u) (M)) < . . . < ben (Nn(M)) < benM. 

(7.6) 

Indeed, JO]) yields ben (Mr W) =J2i<j hen Pi M and benM = ^ P ben P M ; 
which implies (|7.6p . since, by (|6.2p . ben p M is non- negative. 

Definition 8. From now on, we shall note 

7rW=7r»( fl (n)), tt = 7r( ff (n)), £ = £(ff(n)) (7.7) 

so that g{n) — Ntt^ and our work is to compute 7r and £. 

Note that tt and £ are coprime and (|6.5|) implies 

ben gr(n) = ben (Mr£) = ben (Mr) + ben (M;). (7.8) 

Lemma 8. Let j fee a positive integer and S\ < S2 fee two elements of Tj satis- 
fying 

£(5 2 N) <£(5!N). (7.9) 
Then, Si is not the j -prefix of g(n) ; in other words, n^' 7^ Si. 

Proof. If Si = tt^\ equation g(n) = Nt:^ may be written gin) = N (^Si—^j £. 
Set M = N ($2-^y) £ = (S 2 /Si)g{n). From ifO]) . ljT3|) and lTO]l . we get 

£{M) = £(S 2 N) + £ (N^j + £(N£) - 2£(N) 

< £(6iN) + I (^^-y) + t(N£) - W(N) = £(g(n)) < n 

which, from (|1.4p . implies M < g(n) and therefore 82 < Si, in contradiction with 
our hypothesis. Note that our hypothesis implies ben (S 2 N) < ben(<$iAT). □ 
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7.2 Computing plain prefixes 

Let us suppose that we know an upper bound B such that (|6.7p holds. Then 
from (|7.6[) and (|6.7p . for every j such that jx,- < \/x 1 , ben (iW^)) < £> holds. 
Let pj 1 be the largest prime less than ^fx\. Then 7r = n^ 1 ^ (g(n)) is an element 
of T h {B), 

But, we are faced to 2 problems: First, for the moment, we do not know 
B. Secondly, for a given value B' , the sets Tj(B') are too large to be computed 
efficiently. 

What we can do is the following. Let B' < B\. We shall construct two 
non-decreasing sequences of sets Uj = Uj(B') and T>j = T>j(B') with Dj C Uj C 
7~j(B') satisfying the following property: T>j contains the j -prefix tt^ of gin), 
provided that beng(n) < S holds. 

These sequences arc defined by the following induction rule. The only ele- 
ment of T is 1. We set U =V = {1}. And, for j > 1, 

• We define Uj = {5p] \ S G Vj-x, 7 > -a Pj and ben (NSpj) < B'}. 

• By lemma [SJ if 8\ € Uj and if there is a ^2 in Uj such that <5i < 62 and 
£(JV6i) > t{N5 2 ), then 5i is not the j-prefix of g{n). The set V 3 is 
from which are removed these <5i's. In other words, T>j will be the pruned 
set of Uj (see Section |2~2"|) . 

For each 5 in Z>j_i, belongs to if 7 > — a Pj and ben(A' r (5pJ) < B' 
which, according to (|6.5p . can be rewritten as 

ben (Npj) < B' - ben (NS). (7.10) 

It results from Lemma [HI that ben(iVpJ) is non-increasing for — a Pj < 7 < 0, 
non-decreasing for 7 > 0, vanishes for 7 = and tends to infinity with 7. 
Therefore the solutions in 7 of (|7.10p form a finite interval containing 0. 

Thanks to (|7.6p . by induction on j, it can be seen that if ben g(n) < B' , the 
7-prcfix 71"") of g(n) belongs to Uj and also to T>j, by Lemma [5J 

We set "D(B') ~ T>j 1 (B') and since it = Hj 1 , V(B') contains the plain prefix 
7T of g(n), provided that beng(n) < B' holds. 

This construction solves our second problem: at each step of the induction, 
the pruning algorithm makes T>j(B') smaller than Uj(B'), and as we progress, 
Vj(B') becomes much smaller than Tj(B'). 

7.3 Computing B, an upper bound for the benefit 

It remains to find an upper bound B such that (|6.7p holds. The key is Proposi- 
tion [21 Every M such that £(M) < n gives an upper bound for beng(n) + n — 
l{g{n)): 

ben (/(n) < ben g(n) + n — l(g(n)) < ben M + n — £(M). 

We choose some B' , a provisional value of B satisfying B' < B\ . Then 
we compute the set V = T>(B'), and by using the prefixes belonging to this set 
we shall construct an integer M to which we apply Proposition [21 

1 In view of l|6.8jl and after some experiments, our choice is B' = p for 2485 < n < 10 10 
while, for greater n's, we take B' = p/2, and for smaller n's, B' = B± — e where e is some 
very small positive number. 
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Let us recall that pk denotes the greatest prime dividing N. To an element 
S G V(B') and to an integer ui, we associate 

!Spk+iPk+2 ■ ■ ■ Pk+u if w > 

5 if w = 

(5/ {pkPk-i ■ ■ ■ Pk+u+i) if w < and p fc+w+ i > ^/aFT. 

From the definition of prefixes, the prime factors of both the numerator and 
the denominator of 6 G V(B') are smaller than y/xi, and thus smaller than the 
primes dividing the numerator or the denominator of 5 U /S. 

First, to each S G T>. let u = uj(S) be the greatest integer such that ((NS^) < 
n (if there is no such u>(5), we just forget this 5). We call 5^ an element of 
V which minimizes ben (N8 { ° ] ) + n- £(N5 ( ° ] ) and set M = NS { ° ] . From the 
construction of M, we have £(M) < n. By Proposition [H inequality (|6.7|) is 
satisfied with B = benM + n - £(M). 

If B < B\ we stop and keep B; otherwise we start again with B instead of 
B' to eventually obtain a better bound. 

For n = 1000064448, the value of p defined by JSHJ is equal to p w 12661.7; 
the table below displays some values of B' / p and the corresponding values of 
Card(P(£?')) and B/p given by the above method. 



B'/p 





0.2 


0.4 


0.6 


0.7 


0.8 


0.9 


1 


1.1 


\V{B>)\ 


1 


11 


34 


76 


109 


139 


165 


194 


224 


B/p 


7.5 


1.15 


1.13 


1.104 


1.098 


1.082 


1.074 


1.055 


1.055 



In this example, if our first choice for B' is 0.6p, we find B = 1.104p. Starting 
again the algorithm with B' = 1.104/9, we get the slightly better value B = 
1.055p. 

The value of B given by this method is reasonable and less than 10% 
more than the best possible one: for n = 1000366, we find B w 436.04 while 
ben(g(n) + n — £(g(n)) w 406.1; for n = 1000064448, these two numbers are 
13361.6 and 13285.7. 

7.4 How many plain prefixes are there? 

Let us denote by B = B{n) the upper bound satisfying (|6.7[) as computed in 
Section 1731 Let us call n the integer in the range £(N)..£(N') — 1 such that 
B(n) is maximal. 

Let us denote by v — v(n) the number of possible plain prefixes as obtained 
by the algorithm described in Section 17.21 Actually, this number v depends on 
B = B(n) and we may think that it is a non-decreasing function on B so that 
the maximal number of prefixes used to compute g(m) for £(N) < m < £(N') 
should be equal to v(n). 

For the powers of 10, the table of FiG.[5] displays n, n, the quotient of 
the maximal benefit B(n) by p, the maximal number of plain prefixes v(n) 
and the exponent log^(n)/logn. Note that replacing logn by logn will not 
change very much this exponent, since with the notation of Definition 01 we 
have \n-n\< £(N') - £{N) < p k+ i < y/nlogn. 

The behaviour of v(n) looks regular and allows to think that 

v{n) =O(n - 3 ). (7.11) 
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v{n) = # of 


exponent = 


n 


n 


B(n)/p 


plain prefixes 


\ogv(n)l logn 


10 6 


10 3 - 11 


0.9289 


14 


0.3820 


10 4 


10 4 - 10 


0.8453 


19 


0.3197 


10 5 


10 5 - 123 


0.8095 


22 


0.2685 


10 6 


10 6 + 366 


0.9186 


51 


0.2846 


10' 


10 7 - 1269 


0.7636 


59 


0.2530 


io 8 


10 8 + 639 


1.180 


85 


0.2412 


10 y 


10 9 + 64448 


1.055 


212 


0.2585 


10 1U 


10 iU + 88835 


0.6884 


252 


0.2401 


i nil 

10 11 


10 11 + 1007566 


0.9278 


657 


0.2561 


10 i2 


10 12 + 2043578 


1.118 


2873 


0.2882 


IO 13 


10 13 + 5276948 


0.8331 


3805 


0.2754 


10 i4 


10 14 + 17212588 


0.6669 


7048 


0.2749 


10 


10 15 - 44672895 


0.6433 


15148 


0.2787 


1U 


1U — 4oy±zy±y 


n c;n77 




v.zi oy 


10 17 


10 17 - 426915678 


0.6001 


72341 


0.2858 


10 18 


10 18 + 385838833 


0.3027 


144807 


0.2867 


10 19 


10 19 - 9639993444 


0.2963 


170151 


0.2753 


1Q 20 


10 20 + 12041967315 


0.3218 


412151 


0.2808 



Figure 5: The number of plain prefixes. 



7.5 For ben (M) small, prime factors of £(M) are large 

If the number B computed as explained in Section 17.31 is greater than B\ our 
algorithm fails. Fortunately, we have not yet found any n > 166 for which that 
bad event occurs. 

Proposition 4. If B is computed as explained in Section \ 7. 3\ (so that (|6.7p 
holds ) and satisfies B < B\ ( where B\ is defined in (|5.3|) ) then, in view of 
(15. 4p , there exists a unique real number t\ such that 

X\ 

2 < x 2 < Jx{ < p = < ii < xi (7.12) 

loga;i 

and 

p\ogt 1 -t 1 =B. (7.13) 
Further, if ben M < B, we have 

1- Ifx 2 <p< h then v p (M) > 1 = v p (N). 

2. If x 2 < p < y/x~[ then v p (M) e {1,2} and v p (N) = 1. 

3. Ify/x x <p<h then v p (M) = v p (N) = 1. 

4- Ifh<P< xi then v p (M) £ {0, 1} and v p (N) = 1. 

5- Ifxi<p then v p (M) £ {0, 1} and v p (N) = 0. 
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Proof. The function f(t) = plogt — £ is increasing on [x 2 , p] and decreasing on 
[p, x{\. From (|4.8| and (|5.3j) we have 

2 

/(p) > /(a*) = ^ ~ ^ logx 2 - .t 2 = x\ - 2x 2 > B 1 > B > = f(x t ) 

which gives the existence and unicity of t\, which belongs to (/?, x\). Now we 
prove points 1,2,3,4,5. 

Let p be a prime number satisfying x 2 < p < t±. If p does not divide M, 
from (|6.3p and (|6.2p we have 

bcnM > bcn p M = plogp - p = f(p) > f(h) = B. 

Since benM < B is supposed to hold, there is a contradiction and 1 is proved. 

Since we have assumed that B < B\ holds, Proposition [3] may be applied. 
Point 2. follows from point 1. and from item 2. of Proposition [3j while point 3. 
follows from point 1. and from item 1. of Proposition [3l Finally, points 4. and 
5. are implied by item 1. of Proposition [3J □ 



Corollary 1. Let us assume that B is such that (|6.7|) and B < B\ hold. Then 
the suffix £ = £(g(n)) defined in Definition^ can be written as 

t = = p ^--- p ^ u > 0, v > (7.14) 

Pj /'...•••/',. 

where ( we recall that pk is the largest prime factor of N ) 

2 < x 2 < V^i < P < *i < Ph < Pn ' ' ' < Pjv <Pk <Pi x <••• < Pi u ■ (7.15) 

7.6 Normalized prefix of g(n) 

Definition 9. Let u and v be as defined in (|7.14[) and lu = u — v. We define 
the normalized suffix a of g(n) by 

1. Ifio > 

Ph ■■■Pt u £ 



Ph ■ ■ ■ PjvPk + l ■ ■ ■ Pk + U Pk+1 ■ ■ ■ Pk+UJ 

2. Lf lu < 0, we set lu' = —lu and 

Ph ■ ■ -PiuPh ■ ■ -Pk-u'+l , 

O- = = £,Pk ■ ..Pk-u'+l- 

Ph-'-Ph 

The normalized prefix II of g(n) is defined by 

( TTPk+lPk+2 ■ ■ ■ Pk+uj «/w>0 

n =fH * (7-16) 
Na ifw<0. 

\Pk ■ ■ -Pk-u'+l 

Proposition 5. Let a be the normalized suffix of g(n). Then 

Q1Q2 ■ ■ Qs 

a = 

qiQ2 ---qs 

where s is a non-negative integer with 
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1. Ifcj>0 then u < s < v and 

ben (MI) = ben {Ntt)+J2 ben (Np k+i ) = ben (Mr)+^( 2 5 fc+i -plog ] 5 fc+i ), 

i=l i=l 

(7.17) 

= ~Y^{Qi-1i) = Pii+- ■ ■+Piu-(Ph+- ■ -+Pjv)-(Pk+l+- ■ -+Pk+oj) > 0. 

i=l 

(7.18) 

,2. If uj < then v < s < u and, with u>' = —to = v — u, we have 

u)'-\ ( N \ 

ben (MI) = ben {Ntt)+ ^ ben = ben (iV7r)+ ^ (plogp^-pk^) 

»=0 \Pk-iJ i=0 

(7.19) 

S 

%) =X]( < 5 <_% ) =P»i+---+P» u ^(fti+---+ftJ+(? 5 fe+---+Pfe-^'+i) ^ °- 
»=i 

(7.20) 

In 6oi/i cases we have also 

V^T < p < h < qi < ■ ■ ■ < q s < p k+u < Qi < ■ ■ ■ < Q s . (7.21) 
Proof. If u > v then ui = u — v > 0, 

Pii---Pi u £ 



Pjl ■ ■ ■ PjvPk + l ■ ■ ■ Pk + UJ Pk+i ■ ■ ■ Pk + U 



(7.22) 



Since the prime factors pi 1 . . . pi u of the numerator are distinct of the prime 
factors pj 1 . . .pj v of the denominator, a can be written after simplification 

a = QlQ2 --- Qs (7.23) 

9192 • ■ -q s 

where v < s < u and, from (|7.15|) . we have 

^fxl < p <h < qi < q 2 < ... < q a < Pk+^u < Qi < Q2 < ■ ■ ■ < Qs 

which is (|7~2"T|) . From $H>]) we get (f77T7|t while (f7TTH|t follows from ([732]) and 
1(725]) . 

Similarly, if u < v holds, ui' = v — u > 0. So, w' < v, and from (|7.15|) . 
Pk-u'+i > Pk-v+i > P n > h', Q7.22jl and (|7.23|) become 

Pji ■ ■ ■ Pj v 9l • • • Qa 

where u < s < v and we have 

s/x{ < p < ti < 91 < ■ ■ ■ < q s < Pk-u,' = Pk+u < Qi < ■ ■ . < Q s - (7-25) 

which is again (|7.2ip . 

By definition, any prime factor of 7r is smaller than y/xi- Therefore, by 
(|7.25p . Pk-u'+i is greater than any prime factor of n, ()6.5p can be applied and 
(f7TT7]) becomes ((719]) while (f77T8|) becomes ([7720]) . □ 
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The value of the parameter to can be computed from the following proposi- 
tion. It is convenient to set S^, — J2t=i Pk+% (for ui > 0) and S u — — J2i=o 1 Pk-i 
(for oj < 0). In both cases, from (|6.4p . we have 

S u = £(NIL) - £(Nw). (7.26) 

Proposition 6. The relative integer lu which determines the normalized prefix 
II of g(n) (cf. (|7.16|) ) satisfies the following inequalities: 

n - £{Ntt) ?-r- < n - £(Ntt) - - - hCU <S w <n- £(Nir) (7.27) 

1 - p/ti 1 - p/ti 

where 7r is the prefix of g(n) and B and t\ satisfy (|6.7|) and (|7.13D . 

Proof. Let us prove Proposition O for lu > 0; the case lu < is similar. From 
(f7723|) . ([7T2~Tj) and ([738]) . Lemma Q] (i) yields 

1 < a < cxp ( ^ ] • (7.28) 



From ([734]) and ([77T5]) . we have 

£(N£) - £(N) = Ph + ...+ Piu - (p h +...+p jv )= £(a) + S u . (7.29) 
So, we get successively 

ben(Mj) = t(N£) - £{N) - plog£ by (16~T1) 

= + ^{Pk +l - plogPfc +4 ) - ploga by 

i=l 

> e(a)+jT(p k+i -plog Pk+i )-^ by (HMD 

i=l 1 

= £{<j) ( 1 - — ] + ben (MI) -ben (Mr) bv (17471) . 



From (|7.18|) . we have £(<r) > 0. Since, from (|7.21|) . p < h holds, the above 
result together with ([7.8)1 . ()6.7p and (|1.3p implies that 



Q < ^ < ben (Mf) - ben (MI) + ben (Mr) _ ben g(n) - ben (MI) 



1 - p/tx 1 - pit 



^-^(NU) n + £(g(n)) < g - ben (MI) _ _ 
l-p/h l-p/h 

Now, from ([BT3j) . and ([7^5]) . we get 

£(ff(n)) = i(Nn£) = £(Nir) + t(N$) - £(N) = £(Ntt) + £(a) + S u . (7.31) 
Further, since 

n - £{Nn) = £{g{n)) - £(Mr) +n- £(g(n)) = £{a) +S u} +n- £{g{n)), (7.32) 
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we get from (j7~30) and (TO]) 

n - £{Ntt) - B - ben { iVn) < < n - l(Nir) (7.33) 
1 - p/h 

and ([737) follows, since ben (ATI) > 0. Note that (|7. 33[) implies 

ben (MI) < B, (7.34) 

□ 

7.7 Computing possible normalized prefixes 

In Section l?~2l we have computed B such that (|6.7p holds and a set V = V(B) 
containing the plain prefix 7r of g(n). By construction, we know that any prime 
factor of 7r £ T> is smaller than ^fx~i and thus, from (|7.12j) . smaller than t\. 

Definition 10. We call possible normalized prefix a positive rational number 
II = II(7f,w) oftheformll = 9p k+1 . ..p k+UJ (withuj > 0) orll = tt / {p k ■ ■ -Pk+^+i) 
(with ui < 0), where tt € ~D(B) is a plain prefix, and satisfying 

Pk+uo+i > h (7.35) 



n - £{Nn) < n - £(Nt}) - - hCU < < n - £(Nn) (7.36) 

1 - p/ti 1 - p/ti 

with = Yn=iPk+i (ifu>0) and S w = ~Ym=q Pk-i (if u < 0). 

Let us denote by Af the set of possible normalized prefixes; Af has been 
defined in such a way that the normalized prefix II of g(n) belongs to Af. Indeed, 
from (|7.16p . II has the suitable form, the plain prefix tt of g{n) belongs to T>(B). 
flZHU) is satisfied by Proposition and ([735) by (f7~2l) . 

Let us observe that, if uj increases by 1, by (|7.21[1 . S u increases by at least t\. 
In practice, 1 — p/t\ is close to 1 and B is much smaller than t\ so that for most 
of the 7? 's there is no solution to (|7.36[) and there are few possible normalized 
prefixes. For n in the range [998001, 1000000], the number of possible normalized 
prefixes is 1 (resp. 2 or 3) for 1439 values (resp. 547 or 94). For instance, for 
n = 998555, the three possible normalized prefixes are 1, 43/41, 11/10. 

Finally, for a reason given in the next section, for every II G AT, we check 
that the following inequality holds: 

Pk+^+i - (n - £(NU)) > v^T. (7.37) 

This inequality seems reasonable, since, from (|7.35[) . we have Pk+u+i > ti with 
ii close to xi, and, from |73B) . n - £{NIl) =n- £(Nn) - S w < B/(l - p/h) 
which is much smaller than x\. We have not found any counterexample to 

(E33). 
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7.8 The heart of the algorithm 

We have now a list TV of possible normalized prefixes containing the normalized 
prefix II of g(n). For II = 11(7?, u>) G J\f let us introduce 

g{% n) = NUG(p k+ ^,n - £(NU)) = N %9l9lll^l ( 7 . 38 ) 

9i92 • ■ ■ <7s 

where G(p k+lAj ,n - £(iVn)) = Ql< ^ 2 ' " ' — is defined by (021) . We shall use 

9i 92 ■ ■ - 9s 

the following proposition to compute g(n). 

Proposition 7. The following formula gives the value ofg(n): 

g(n) = max g{fl 7 n) = max NUG(p k+UJ , n - £{Nn)). (7.39) 
neW neW 

Proof. Note that (|1.13|) and (| 1 . 14[) imply either s = or the smallest prime 
factor q s of G(p k + UJ ,n — £(NIi)) satisfies pt+u+i ~ Qs < n — l(NH) which, 
from (|7.37|) . implies q s > ^fxi and thus, the prime factors of ir and those of 
G(pk+u,n — £(NH)) are distinct. Therefore, for any II = H(tt,uj) g TV with 
u > 0, we get from ((733), JB3J and (TlTT5|) 

«(ff(3,n)) = 1(JVtt) + £ f 7V Pfc+1 " ' ' Pfc+ " Ql ' ' ' Qs ) - £(N) 

V qi- ■ - q s J 

UJ s 

= £(m) + Y^Pk+i + ^2(Qi-qi) 

i=l i=l 

= «(Wn)+*(G(p fc + w ,n-*(iVn))) 
< £(/ATI) +n- ^(-/Vft) = n . 

Inequality £{g(Ii, n)) < n can be proved similarly in the case w < 0. 

Since £(#(±I,n)) < n holds, l|1.4p implies for all ft G Af 

ff(n,n) < flr(n). (7.40) 

From (|7.16p . we get g(n) = ATIer where II is the normalized prefix of g(n). Now, 
if u > 0, from 1(715)) . 1(731)1 . ((TIB)) and ifO)) . we have 

i=i i=i 

= £(g{nj) - l(NH) <n- £{NU) (7.41) 

(£(o-) <n - (.(Nil) still holds for u < 0). Therefore, in view of |7!2l)l and of 
Definition 1)1. 12|) of function G, we have 

g(n) = NTla < NTlG(p k+u ,n - £(NIl)) = g(Il, n). (7.42) 

Since n G TV, 1(712")) and (fflO)) prove 1(735)1 ■ □ 
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7.9 The fight of normalized prefixes 

Let LTi and LT2 two normalized prefixes. By using Inequalities (|8.4p below, it is 
sometimes possible to eliminate LTi or LT2. 

Indeed, from (|8.4|) . we deduce a lower and an upper bound for 17(11, n) (de- 
fined in (17381) ): 

g'(U,n) < g{%n) <g"{%n). 

If, for instance, g"(JI\,n) < <?'(Il2,n) holds, then clearly ITi cannot compete in 
(|7.39p to be the maximum. 

By this simple trick, it is possible to shorten the list TV of normalized prefixes. 
For instance, for n = 10 15 , the number of normalized prefixes is reduced from 9 
to 1, while, for n = 10 15 + 123850000, it is reduced from 37 to 2. 



8 A first way to compute G(pk, m) 
8.1 Function G 

In this section, we study the function G introduced in (jTTTSJ) . First, for k > 3 
and < m < Pk+i — 3, we consider the set 

9(pk,m) = \f= QiQ*~-Q> . e(F) = Y j (Q i -q i )<m, a > o) (8.1) 

where the primes Qi, Q 2 , ■ ■ ■ , Q s , qi, qi, satisfy (|1.13|) . 

The parameter s = s(F) in (|8.ip is called the number of factors of the 
fraction F. If s = 0, we set F = 1 and £(F) = so that G(pk,m) contains 1 
and is never empty. The definition (I1.12[) can be rewritten as 

G(p k ,m)= max F. (8.2) 

Obviously, G(pk,m) is non-decreasing on m and G(pk,2m + 1) = G(pk,2m). 
Note that the maximum in (|8.2[) is unique (from the unicity of the standard 
factorization into primes). It follows from (|1.13|) that, if < m < Pk+i — Pk, 
the set Q{j>k-,fn) contains only 1, and therefore, 

< m < p k+ i - p k G(p k ,m) = l. (8.3) 

Proposition 8. 1. Let q be the smallest prime satisfying q > pu+i — 771. The 
following inequality holds 

Pk+i . nl \ Pk+l (o * \ 

<G{pk,m)< . (8.4) 

q pk+i - m 

Note that if q — Pk+i — m is prime, then (|8.4|) yields the exact value of 
G{pk,m). 

2. Now, let F = ^ ^ ^ be any element of Q(pk,m); we have 
qi 92 ■ ■ ■ q s 

G( Pk ,m)>F>l + ^l- (8.5) 
Pk 
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Proof. The lower bound in (|8.4p is obvious. Let us prove the upper bound. If 
< m < Pk+i —Pk, the upper bound of (|8.4|) follows by (|8.3|1 . If m > p/c+i — Pfc, 
— — g GiPkifn) and thus G(pk,m) > + > 1. Moreover, with the notation 

(|8.ip . if G(pk, m) = F = < ^ 1 ^ 2 — , we have s > 1 and Lemma [T] (ii) implies 

9i 92 • • • q s 

G{p k ,m) < Qs < < (8.6) 

Q s -i[F) Q s -m pk+i-m 

where the last inequality follows from (|1.13p and the decrease of 1 1— > t/(t — m). 

Let us now prove (|8.5p . This inequality holds if t(F) = (i.e., F = 1 and 
s = 0). If s > 0, from |[L~T5]) . we get 

Qi -, i Qi 9i ^ -t i Qi — 9i • 1 o 

— = H >H , i = l,2, ...,s 

qi qt Pk 

and 

F=n^>nfi+^)>i+^^=i+^- 

□ 



8.2 Function i7 

Let M < Pk+i — 3; we want to calculate G(pk,m) for < m < M. Let us 
introduce a family of consecutive primes Po < Pi < ■ ■ ■ < Pk = Pk < Pk+i < 
. . . < Pr < Pr+i (so that Pi = Pk+i-K for < i < R + 1) with the properties 

P R+1 -P K >M, R>K+1, P k+1 -P >M, P x >3. (8.7) 

It follows from (|8.ip and (|1.13p that the prime factors Q%, . . . , Q s , q%, . . . , q s of 
any element of G(pk, m) = G(Pk, tu) should satisfy 

Pi < q s < ■ ■ ■ < qi < Pk = Pk < Pk+i < Qi < ■ ■ ■ < Qs < Pr- (8.8) 

Of course, in ()8.7[) we may choose Pr (resp. Pi) as small (resp. large) as possible, 
but it is not an obligation. 

Let us denote by Q\, Q' 2 , . . . , Q' R _ K _ S the primes among Pk+i, ■ • ■ , Pr which 
are different of Qi, . . . , Q s ', we have 

Q[ + Q' 2 + ... + Q' R _ K _ S = P K +i + ... + Pr-(Qi + ... + Q s ) (8.9) 

and (18.21) becomes 



ni-n \ Pk+xPk+2 ■ ■ ■ Pr Pk+\Pk+2 ■ ■ ■ Pr , S1(1 > 

G(P K ,m) =max— — = . , , ; r (8.10) 

Qi ■ ■ ■ Qr-K-sVi ---qs min( 9l . . . q R _ K ) 

where the minimum is taken over all the subsets {q^, q' 2 , . . . , q' R _ K } of R - K 
elements of {Pi, . . . , Pr} satisfying from (|1.14p and 



q[+q' 2 + ... + q R _ K = Q[ + Q' 2 + . . . + Q' R _ K _ a + qi + q 2 + ■ ■ ■ 

S 

= P K +i + Pk+2 + ■ ■ ■ + Pr - V (Qi - q t ) 



i=i 



> P k +i+Pk+2 + ■■■ + PR-m. (8.11) 
(Note that, from (|87T|) . R - K > 1 holds). 
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Definition 11. For l<r<R,l<j< mm(r,R — K) < R and m > 0, we 

define 

H(j,P r ;m)=xmn(q' 1 q' 2 ...q' :j ) (8.12) 
where the minimum is taken over the j-uples of primes (g^, q' 2 , ■ ■ ■ , q'j) satisfying 

P\ < q[ < q' 2 < ■ ■ ■ < q'j < Pr (8.13) 



q[ + q' 2 + . . . + q'j > P K +i + P K +2 + ■■■ + P K +j - m. (8.14) 
// there is no (q[, q 2 , ■ ■ ■ , q'j) such that (|8 . 1 3[) and (|8.14p hold, we set 

H(j,P r ;m) = +oo. (8.15) 

By the unicity of the standard factorization into primes, the minimum in 
(PSH) is unique and ([Q]]) and <|HTT^|) yield 

G(p k ,m) = G(P K ,m) = ^f K ; 2 p "^ (8.16) 

H [R — K, Pr ] m) 

For j = R — K and r = R, the j-uple q[ , q' 2 , . . . , q'j defined by q\ = PK+i satisfies 
^Wf and ([HH]) for all m > 0; so, H(R—K, P R ; m) is at most P k +iPr+2 ■■■Pr 
and is finite. 

8.3 A combinatorial algorithm to compute H and G 

Definition 12. For every integers (r,j), 1 < r < R and 1 < j < R — K, we 

define 

m,j(P r ) = 

f P K +i + Pk+2 + ■■■ + P K +j - (P r + P r -i + ■ ■ ■ + P r -j+i ) ifj<r 
1 +oo if j > r. 

(8.17) 

Remark: If j > r+1, (|8.13p cannot be satisfied and, from (|8.15|) . H(j, P r ; m) = 
+oo for all m > 0. If j < r, from (|8.14|) . it follows that, if m > mj(P r ), 
H(J, P r ; m) < P r P r _i . . . Pr-j+i while, by |8^l5]l . if m < m 3 (P r ), H(J, P r ;m) = 
+oo. So that, in all cases, if to < m.j(P r ) 1 H(j, P r \ m) = +oo. 

Note that, for j fixed, nij(P r ) is non-increasing on r since, for j < r, 

p p sn -f ^ v i (8 ' 18) 
P r — P r -j > if 1 < J < r — 1, 

and, for j > r + 1, mj(P r ^i) and mj{P r ) are both +oo. On the other hand, if 
j < min(r, R — K) for every to such that 

to > Mj{P r ) = P K +i + Pk+2 + ■■■ + P K + 3 - (Pi + P 2 + ■ ■ ■ + Pj), 
H(j, P r ; to) is equal to P\P 2 ■ ■ ■ Pj- 
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Proposition 9. For j = 1, from (jHTT^jl . (j5TT3"|) and (fgH|) . we Ziaue 



H(l,P r ;m) 



(8.19) 



fP : ifm>M 1 (P r )=P K+ i-Pi 

Pi ifl<i<r and P K +i - Pi < m < Pk+i - Pi-i 

P r ifmi(P r ) = P K+ i - P r < m < P K +i - P r -i 
oo ifm < mi(P r ) = Pk+i - P r - 

Further, we have the induction formula: 

H(j, P r ;m) = min (H(j, P r _ i; m), P r H(j - 1, P r _ i; m - P* +j + P r )) . (8.20) 

Proof. The calculation of H(l, P r \ to) is easy. Let us show the induction formula 
([cT2"0)) . Either P r does not divide H(j,P r ;m) and H(j, P r ;m) = H(j , P r _i;m) 
or P r = is the greatest prime factor of H(j, P r ;m) = q\q! 2 . . . q j and from 
(EH, we get q[ + ...+ > P A -+i + . . . + Pjr+j-i - (m - P A - +J + P r ). □ 

Note that if to > rrij(P r ), m — Pk+j + P r > mj_i(P r _i) since m 3 (P r ) = 
TOj-i (Pr-i) + P/f+j ^ Pr so that if (j, P r ; to) and if (j — 1, P r _i; m — Pic+j + Pr) 
are simultaneously finite or infinite. (|8.18jl implies that mj(P r ) and m,j{P r -\) 
are both infinite or mj(P r -i) > mj{P r ). For m,j(P r ) < to < ?7ij(P r _i), (|8.20p 
reduces to 

if(j, P r ; to) = P rJ ff(j - 1, P r _i;m - P A - +J + P r ) (8.21) 

while, for to > mj(P r _i), the three values of the function if in (|8.20p arc finite. 
From (|8.19| . we may remark that, if we set 



if (0, P r ; m) = 1 for all r > 1 and m > 0, 



1.22) 



the induction formula (|8.20j) still holds for j = 1 . 

In view of ([8TT5]) . for 1 < r < R, 1 < j < min(r, R - K) and mj(P r ) < 
to < M, we calculate if(j, P r ; to) by induction, using for that (|8.22p . (|8.20p and 
(f8~2"Tj) . If if + 2 < r < R, it is useless to calculate H(j, P r ; m) for j < r - A'. 

Finally, after getting the value of H(R — K, Pr;to) for tor_a-(Pr) = < 
to < M, we compute G(pk,m) by (|8. 16|) . 



8.4 Bounding the largest prime 

It turns out that the largest prime used in the computation of G(pfc,m) for 
< m < M is much smaller than Pr defined in (18. 7|) . For instance, for p,t = 
P A - = 150989 and A/ = 5000, P defined by (jE2]) is at least equal to K + 425 
while only the primes up to Pk+5 = Pk+5 = 151027 are used. 

So, the idea is to replace R by a smaller number R, K + 1 < R < P, and to 
calculate by induction if (P — K, Pg; to) instead of H(R — K, Pr; to). We get 

the fraction F = K t} — } S22l R which satisfies F < G(pk, to). Now we have 

H(R- K,Pg,m) 
the following lemma. 
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Lemma 9. Let F be a real number satisfying 1 < F < G(pk, to) 



9i?2 •• -q s 

Then, the largest prime factor Q s of the numerator of G{pk, rn) is bounded above 
by 

( mF \ 

Qs < min U fc + to, y—j J • (8.23) 

Proof. Using Lemma 1 and ()1.15p . we write 

c ^ ri, \ Q1Q2 • • • Qs . Qs ^ Qs 

F < G{p k ,m) = < — — — — < 



qiq 2 ...q s Q s ~ £(G(p k , m)) Q s -m 

mF 

: ^ yields Q s < . On the other hand, Inequality (|1.13|) together with 

F — 1 

(|1.14p implies Q s — Pk < Qs — 9s < m which completes the proof of (|8.23|) . □ 

IfF = P k+iPk+2 ■ ■ ■ Pr > x . f > / + mF_ \ ^ foUows 
H(R-K,P n ;m) R \ F-lJ 

from Lemma IHl that G(pk,m) = F. If not, we start again by choosing a new 

( m F \ 

value of Pg greater than min Pk + m, — . Actually, Inequality (18.231) 

\ F - l) 

gives a reasonably good upper bound for Q s . In the program, our first choice is 
R = K + 10. 

8.5 Conclusion 

The running time of the algorithm described in sections 18 . 31 and IS^O to calculate 
G(p,m) for m < M grows about quadratically in A/, so, it is rather slow when 
M is large. 

For instance, the computation of (?(10 15 — 741281) leads to the evaluation of 
G(p, 688930) for p = 192678883, and this is not doable by the above combina- 
torial algorithm. 

In the next section, we present a faster algorithm to compute G(pk, rn) when 
to is large, but which does not work for small ?«'s so that the two algorithms 
are complementary. 

9 Computation of G(pk,m) for m large 

The algorithm described in this section starts from the following two facts: 

• if G(pk, rn) = ^1 anc ] m j s large, the least prime factor q s of the 

denominator is close to pk+i—rn while all the other primes Qi, . . . , Q s , qi, . . . 

q s -i are close to pk- More precisely, G(pk,m) is equal to +1 G(pk+i, d) 

Qs 

where d~m — Pk+i + <Zs is small. 

Note that when m is small G{pk, rn) is not always equal to 

^^G(p k +i,m - Pk+i +q s )- For instance, G(103, 22) = * while 
q s 97 x 101 
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• In (|8.5p . we have seen that £(G(p, to)) = to implies G(p, to) > 1 + and 
it turns out that this last inequality seems to hold for to large enough. 



9.1 A second way to compute G(pk,m) 

We want to compute G(pk,m) for a large to. The following proposition says 
that if, for some small 5, pk — m + 5 is prime and such that G(pk+i, 5) is not 
too small, then the computation of G{pk,rn) is reduced to the computation of 
G{pk+\, to') for few small values of to' . 

Proposition 10. We want to compute G(pk,m) as defined in (|1.12|) or (|8.2p 

wii/i pfc odd and Pk+i ~ Pk < to < Pk+i — 3. VFe assume that we know some 
even non-negative integer 8 satisfying 

Pk+i + 5 — m is prime, (9-1) 
G(p k+1 ,6)>l + — (9- 2 ) 

and 

(9-3) 

If S = 0, we know from Proposition^ that G(pk,m) = — P k+1 jj § > 0, we 

Pk+l - TO 

have 

G(p k ,m)= max — — G(p k +i,m - p k +i + q), (9.4) 

q prime q 
Pfc + l-m < q < q 

where q is defined by 

_ p k +iPk+2(Pk+i ~m + 8) 35 

9 = 7 — ?w tTavT <Pk+2-m+— ■ (9.5) 

CPfe+i + d ){Pk+i - 3(3/2) 2 

Before proving Proposition |TD] in Section 19. 3[ we shall first think to the 
possibility of applying it to compute G(pk, m). 



9.2 Large differences between consecutive primes 

For x > 3, let us define 

A(x) = max(pj —pj_\). (9-6) 

Pj<X 



Below, we give some values of A(x): 



X 


10 2 


10 3 


10 4 


10 b 


10 b 


10 v 


10 y 


10 y 


10 i0 


10 11 


10 iy 


A(x) 


8 


20 


36 


72 


114 


154 


220 


282 


354 


464 


540 


(logx) 2 


21 


48 


85 


133 


191 


260 


339 


429 


530 


642 


763 



A table of A(x) up to 4- 10 12 calculated by D. Shanks, L.J. Lander, T.R. Parkin 
and R. Brent can be found in [26], p. 85. There is a longer table (up to 8 • 10 16 ) 
on the web site [TB]. H. Cramer conjectured in 4] that lim^^oo fa^yi = 1- For 
x < 8 • 10 16 , A(x) < 0.93(loga;) 2 holds. 
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Let us set A = A(pk+i); let us denote by Si = Si(pk) the smallest even 
integer such that i5i > A and 

G(p k+1 ,d)>l + —, d = S 1 -A + 2,S 1 -A + 4,...,S 1 . (9.7) 
Pk+i 



By using the combinatorial algorithm described in 18.31 we have computed that 
for all primes p k < 3 • 10 8 , we have 5i(p k ) < 900 = <5i (252314747) and 

5i( Pk ) < 2.55(logp fe ) 2 . (9.8) 

To compute the suffix of g(n) for n < 10 15 , we do not have to deal with larger 
values of Pk- However, for larger p k s, we conjecture that 5i(pk) exists and is 
not too large. 

Lemma 10. Let p k satisfy 5 < pk < 3 • 10 8 , m be an even integer such that 
Pk+i - Pk <m< p k+ i - 3, and Si = Si(p k ) defined by (j!JT)) . If m > |5i(p fe ), 
then there exists an even non-negative integer 

5 = S(p k ,m) < Stipk) < 2.55(log Pfc ) 2 (9.9) 

such that (|9.1I) . ()9.2jl and (|9.3j) hold. Therefore, Proposition \l(A can be applied 
to compute G(pk,m). 

Proof. Let us set a = Pk+i + <5i(j>fc) — m. We have 

7 7 

a = Pk+i + Si(pk) - m < pk+i - -Si (pk) < Pk+i - - A < p k+1 . 

Since 5\ > A and m < pk+i — 3, a > A + 3 holds. From the definition 
of A = A(pk + i), there exists an even number b, < b < A — 2 such that 
a — b = Pk+i — m + (Si — b) is prime. From the definition of S\{jpk), we know 
that G(p k+1 ,S 1 -b) > 1+ Therefore, S = S 1 -b satisfies (^TJ) , (JOJ), ^1 

and < S < S^pk). The last upper bound of (HU) follows from □ 

9.3 Proof of Proposition [T01 
A polynomial equation of degree 2 

Lemma 11. Let us consider real numbers Ti,T2,S satisfying 

< Ti < T 2 (9.10) 

and 

IT 

(5 = or 5>T 2 -T 1 ) and S< — -■ (9.11) 



Note that (f^TTTT]) and ((9~TT|) imply 



Let m be a parameter satisfying 



T 1+ S<-^^-- (9.12) 



2 

We set 



9S 

< — < m < T x . (9.13) 



E(X) = X 2 -(T 1 +T 2 -m)X+?^±-j r -^- (9.14) 

±1+0 
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1. The equation E(X) = has two roots Xy and X2 satisfying 

< Xy < Tl+T ^~ - <X 2 <T 2 -8. (9.15) 



2. For Ty,T 2 and 8 fixed and m in the range (|9.13|) . X 2 is a non- decreasing 
function of m. 



3. We have 



T 1 -^< T -±±^l-^<X 2 <T 2 -8. 
2 3 2 ~ 



4- Let Yy and Y 2 be two positive real numbers satisfying 
Yy < Y 2 , Yy+Y 2 =Ty+T 2 -m and > 



Y{Y 2 ~ Ty + 8-m 



(9.16) 



(9.17) 



We have 



Y 2 > X 2 > T n - ? 



Q SJ 

Y x < X x < T 2 -m+ —. (9.18) 



Proof. 1. The discriminant D of (|9. 14)) can be written as 

n9 TyT 2 (Ty + 8-m) 

d = r 1 + r 2 -m 2 -4 21 1 

ly + 
4:8 m 2 T ? 



{m + T 2 -Tyf 



m[m + T 2 -Ty) 2 (Ty + 5) 



(9.19) 



since, from (|9.10p and ()9.13p . m + T 2 — Ty docs not vanish. If 8 = 0, the 

T 2 

above bracket is 1 while if 8 > T 2 — Ty > 0, the fractions — and 



Ty +8 



m + T 2 — Ty 



are at most 1, so that in both cases (|9.19|) yields 



D > {m + T 2 -Tyf 



48 



(9.20) 



Therefore, from and (r9~TUl) . D > ( m + T2 — I±L > holds. 

9 

The sum Xy +X 2 of the two roots is Ty +T 2 — m which explains the second 
and the third inequality of (|9.15p . Further, since Ty < T 2 and m > 28, 

Tl + ^ 2 ~ m < T 2 - 8 holds. By (EH, and (l9~T2l . 



S(T 2 - <5) = (Ti + <5 - to) 



(T 2 - 6) > 



which proves the last inequality of (|9.15D . 

Remark: If 8 = 0, the roots of (pHi)) are Xy = Ty - m and X 2 = T 2 . If 
8 = T 2 — Ty , they are Xi = T 2 — to and X 2 = Ti . 
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2. By (19.141) . X 2 is implicitcly defined in terms of m and, through (|9.12|) . we 

have 

dX 2 _ -$% _ Zfe-Xt > T 2 -S-X 2 



dm |f 2X 2 -(T 1 +T 2 -m) - 2X 2 -(T 1 +T 2 -m) 
which is non- negative from (|9.15p . 

3. For m=f, dOfl yields y/D > m + T ^ = f + ^II and 

„ Ti + T 2 - m + VD ^ Tx + 2T 2 36 36 

A9 = > — > 1 1 — 

2 3 2 ~ 2 

Further, for m > the upper bound in (|9.16|) follows from (ii). 

4. Conditions IPTt imply Bfyj = £(Y 2 ) = -y 1 y a + riT2 ^ 1 + ^ ~ m ^ > 

Ti + 

so that Yi < Xi and Y 2 > X 2 ; (f^TTH|) follows from (|9~TC|) and from X x = 

Ti + T 2 - m - X 2 . 

□ 

Structure of the fraction G(pk,m) 

Lemma 12. Let k and m be integers such that k > 3 and Pk+i — Pk < m < 
Pfc+i — 3. VFe write 

G(p fc ,m)=F = Qlg2 --- Qs (9.21) 
9192 • ■ - 9s 

m'tt s > 1 cmd Qi, . . . , Q s , qi, . . . , q s primes satisfying 

3 < q s < q s -i < ■■■ <qi <Pk< Pk+i < Qx < ■ ■ ■ Qs-1 < Qs, (9.22) 



Pk+ 



i-Pk< ((F) = ^(Qi - qi ) <m< Pk+1 - 3 < Pk+1 (9.23) 



(=1 



and we assume that there exists an integer 6 such that 
2m 

0<6<—, and (5 = or 6 > Pk +2 - Pk+i) (9.24) 

and 

F> Pk+1 + S ■ (9.25) 
Pk+i -m + 5 

We apply Lemma \Tl\ with T\ = Pk+i and T 2 = Pk+2, 6 and m, and we denote 
by X\ and X 2 the two roots of equation (|9.14[) . E(X) = 0. Then we have 

1- Q s <Pk+i + 6, 



def 

2. for s > 2 and 1 < i < s - 1, K = Qi - qi < Pk+2 - X 2 , 



3. for s > 2 and 1 < j ' < s — 1, Aj === < Pk+2 — X; 

1=1 
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Moreover, if we write F = UV with 

TT Q\Qi ■ ■ ■ Qs-iQs , t t Pk+i ,„ „m 

U = and V = , (9.26) 

qiq 2 ■ ■ ■ q. s -iPk+i q s 

we have, for s > 1 

3(5 

4. £(U) = A s _i + Q s - p k+1 < p k+2 - X 2 < Pk+2 - Pk+i + y 



an 



d 



K , , ~ Pfe+iPfe+2(pfc+i - m + <5) 

(Pk+i + S)(p k+1 - 36/2) 

Proof. 1. First, we observe that (|9.22p implies 

Qi>Pk+i>Pk+i, l<i<s. (9.27) 

Lemma Q] and (|9.23p yield respectively F < g %Tf1 an< ^ ^(^) — TO ' so 
that, together with (|9.25j) . we get 

Pk+i +5 < f < ® s < 



Pk+i +6 -m Q s ~ 1{F) Q s - m 

which, with the decrease of 1 1— > rzjr, gives Q s < p/c+i + 5. 

2. From the definition of Xi and (|9.22[) . A^ is positive and increasing on i, and 
it suffices to show A s _i < pk+2 — X 2 . Wc write F = F\F 2 with F\ = 

and F 2 = n^ s _i f^- From <P^1> and (E22J), we have 

Pk+i > m> m - X s -i > £(F) - A s _i = A x + . . . + A s _ 2 + X s > X s > A s _i 

which implies 

Pk+2 - Xs-i > Pk+i - Xs-i > Pk+i - (m - A s _i). (9.28) 

Further, Lemma [TJ (|9.23|) and the increase of t t— > ^ s _ t , (|9.27p and the 

decrease of t ^ t -(m-\ s ~Tj > im P l y 

T-l . Qs QS 



h-£(F 2 ) Q s - (£(F) - A s _i) 

< ^ < (9.29) 

Q s - (m - A g _i) pfe+i - (m - A s _i) 

If s > 3 or Qi > Pfc+2, (|9.22p implies > which yields i*i = 



< 



Pk + 2 



<2s-l-A s _l — Pfc + 2-A 



- so that, from ((!L2"5|) and ([9~2T?|) , we get 



^+i + <5 < F = F,F 2 < Pk+2 ^ r (9.30) 



Pk+i + 6-m p k+2 - A s _i pk+i - (m - A s _i) 

Let us set Y" 2 =Pfc +2 - A s _i, Yi = p fc+ i -(m-A s _i); from (|9.28j) . F 2 > *i 
holds and, in view of (|9.30p . we may apply Lemma fill Point 4. to get 
Y 2 = Pk+ 2 — A s _i > X 2 which implies 2.. 
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If a = 2 and Q x = p k+u F = = ^ Qa .g. a) . From (D) we 

have Q 2 > Pk+2 and F < Pk +^ pt;+2 ^(Q 2 2 _ g2 ) ■ Here we set F 2 = <7i and 
*i = Pk+2 - (Q2 - q%) =q 2 - (Q2 - Pk+2); by $M and $7%$,, we get 

Y 2 = qi > q 2 > Yi = q 2 - (Q 2 - Pk+2) = Pk+2 - A 2 

2 



> Pk+2 - J~] Xi = Pk+2 - £(F) > Pk+ 2 - m > 0; 



i=i 

we may still apply Lemma [TT] Point 4. to get Y 2 = qi — p k +i — ^i > X 2 . 
which implies 2.. 

3. This time, we write F = F X F 2 with F x = n*=i ff and ^ 2 = IIi=j+i ff 
so that £(^i) = Aj and ^(F 2 ) = £(F) - Aj < m — Aj. For 2 < j < s - 1, 
from (|5^5]l . Lemma U ijOTj) , and we get 



p fc+ i + <5-m- " Qj-l{F x ) Q s -l(F 2 ) 

Pk+2 Pk+1 



< 



Pk+2 - Aj pk+i - (m - Aj) 



Therefore, we apply Lemma [TT] Point 4., but we do not know whether 
Pk+2 ~ Aj is greater than pk+i — {rn — Aj), so that, cither 

Pk+2 - A 3 > X 2 (9.31) 

or 

Pk+2-AjKXx. (9.32) 

For j — 1, as Ai = Ai, (|9.3ip holds, from 2.. Since Aj is increasing on j, 
if (|9.31|) holds for some j = jo , it also holds for j < jo . If (|9.3ip holds for 
j = s — 1, 3. is proved; so, let us assume that the greatest value jo for 
which (|9.3ip holds satisfies 1 < jo < s — 1; we should have 

p k+ 2 - A jo > X 2 and p k+2 - A J0+1 < X x . (9.33) 

From 2., (|9.33p and because X x , X2 arc solutions of (|9.14| . wc should get 

Pk+2 - X 2 > Aj 0+1 = Aj 0+ i - Aj > X 2 - Xi = 2X 2 +m- p k +i - Pk+2 

which, would imply m < 2pk+2 + Pk+i — 3A" 2 and, through the second 
inequality of (|9. 16|) . m < in contradiction with (|9.24| . Therefore, 
jo > s — 1 and 3. is proved. 

4. If s = 1 wc have to show l(U) = Q\— Pk+i < Pk+ 2 — X 2 which is true since, 
from 1., Qi - pk+i < S and from (|9.16p . with T 2 = Pk+2, $ < Pk+2 ~ X 2 . 

So, we assume s > 2. If Qi = Pk+i, U simplifies itself; and, in all cases, 
from (|9.22p , the prime factors of the numerator of U are at least Pk+2 and 
those of the denominator are at most p k +i- So, we may apply Lemma [1] 
which, with (|9.27p and the decrease of t h- ► t/(t — £(U), yields 

U < 9l < **±= V = (9 34) 

~ Q s - £{U) ~ pk+2 - £(U) ' Pk+i - i{V) 1 ' 
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It follows from §£25§ that £{U) + 1{V) = £(F) < m and, from we 
get 

Pk+1 +S < =U y < Pk+2Pk+l 



Pk+1 + S-m~ ~ { P k+2 - t(U))(pk+i - (m - £(U))) 

Applying Lemma [TT] Point 4. with (Yi,Y 2 ) = (Pk+2 — £(U),p k +i — (to — 
t{U))) yields 

p k+2 -t(U)> X 2 or p k+ 2 ~ £(U) < Xl (9.35) 

But, from 1. and 3., we have £{U) = A s _i + Q s — Pk+i < Pk+2 — X 2 + 5 
which, together with (Xi,X 2 ) solutions of (|9.14|) . the second inequality in 
([9~To) and l[9~33[) . give 

Xi + £(U) - Pk+2 < X\ — X 2 + S = 5 + Pk+i + Pk+2 — to — 2X 2 

2 

< S + p k+1 + p k+2 to -(pk+i + 2p fc+2 ) + 35 

AX 1 "Pk+1 - Pk+2 „ 

= 4o H to < 0. 

o 

Therefore, — ^"(l 7 ) < -Xi does not hold, and, from (|9.35[) . we have 

Pk+2 — £{U) > X 2 which shows the first inequality in 4.. The second 
inequality comes from (|9.16p . 

5. From (f^22|) and [pQ3]) . wc have £{V) = p k+1 - q s <Q s -q s < i(F) < m 
which proves the lower bound of 5.. 

If s = 1 and Qi = p fc+1 , [7 = 1 and F = V so that, from ([935)) , 

- Pfc+1 < Pfc+l(Pfc+l ~ TO + 8) Pk+lPk+2{Pk+l - TO + (5) 

qs ~ F ~ Pk +i+S ~ q ~ {pk+i+S)(p k +i-'3S/2)' 
If s > 2 or Qx > p fe+2 , ([Oil) holds and gives with (pOS) and 4. 

_ _ Pk+lU Pk+lPk+2(Pk+l - TO + 8) _ 

K F ~ (pk+i + S)( Pk +2 - £(U)) ~ Q - 

□ 

Proof of Proposition [lOl 

Let us assume S > 0. (|9.2|) and ()8.3[) imply 

^ > Pfe+2 -Pk+i- (9.36) 

First, we prove the upper bound (|9.5p . We have to show that the quantity below 
is positive: 

(Pk+2 - TO + 5)(p fc+ l + 8) [pk+1 - y) - Pk+lPk+2(Pk+l - TO + 8). 

But this quantity is equal to 

(Pk+2 - p k +i)({pk+i - 8){m- y) + S(m - 3S)J 

5( 98\ 3(5 2 / 38 \ 

+ — TO H TO 

m 2 V 2 / 4 V 2 / 
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which is clearly positive since, from (|9.3p . p k +i > m > ^ holds and (|9.5p is 
proved. 

Let q be a prime satisfying Pk+i — m < q < q. In view of proving (|9.4D . let 
us show that 

?t}±G(p k+1 ,m - p k+1 +q)< G{ Pkl m) (9.37) 

q 

holds. Let q' be any prime dividing the denominator of G(p k +\ 1 m — Pk+i + <?); 
we should have Pk+2 — q 1 < m — Pk+i + q i.e., q 1 > Pk+i + Pk+2 — m — q which 
yields from (HH), (pHie]) and d9T3j) 

q' - q > Pk+i + Pk+2 -m-2q> p k+1 + p k+2 -m-2q 

( 3A 

> p k +i + Pk+2 — m — 2 I pk+2 - m + — J = pk+i - Pk+2 +m-35 

> Pk+i - (S + p k +i) + m — 35 = m — AS > 0. 

Therefore, q' ^ q, and after a possible simplification by p k +i, ^^-G(p k +\, m — 
p k+1 +q) e Q(pk,m) (defined in (|HTT|) ). which, from implies ([Q7| . 

From (|9.36[) and (|9.3p . we have < 25 < m, and the prime p = Pk+i + 
5 — m satisfies p < Pk+2 — 8, and thus is smaller than any prime factor of the 
denominator of G(j>k+i,6). Therefore, after possibly simplifying by p k +i, the 
fraction $ = ^ ±L G(p k +\ 1 5) belongs to G(p k ,m) and we have from (|8.2j) and 

urn 

G(p kl m) > $ > — 1 



p k+1 + d-m\ p k +i J p k+1 +d-m 

So, hypotheses (|9.24p and (|9.25| being fullfilled, we may apply Lemma IT2| (v) 
which, under the notation (|9.26[) . asserts that 

G( Pk ,m) = UV = C/^ii (9.38) 

q.i 

with q s G [pk+\ - m,q] and £(U) + £(V) = l{G(jp k ,m)) which, from (fl~T5|) . 
implies ^(t/) < m — ^(V) = m — p k +i + q s - After a possible simplification by 
Pfc+i, L/ belongs to £(p fc+ i, £([/)) C G {p k+1 ,m - p k+1 + q s ). So, from (JO}, 
^ < G(p k+ll m - p k+ i + q s ), and (]9.38[) gives 

G(p k ,m) < ^^-G(p k+1 ,m - p k+1 + q s ) 

qs 

which, with (|9.37[) . completes the proof of (|9.4[) and of Proposition [TUJ □ 



10 Some results 

With the maple program available on the web-site of J.-L. Nicolas, the fac- 
torization of g(n) has been computed for some values of n. The results for 
n = 10 6 , 10 9 , 10 12 , 10 15 are displayed in Fig. 6. For primes qi < q-2 let us de- 
note by [91 — 92] the product ri gi <p< 92 P- The bold factors in the values of g(n) 
are the factors of the plain prefix ir of g(n), defined in ([$]). 

On a 3GHz Pentium 4, the time of computation of g(n) is about 0.02 second 
for an integer n of 6 decimal digits and 10 seconds for 15 digits. 



39 



n = 10 6 , N = 2 9 3 6 5 4 7 3 [ll-41] 2 [43-3923] 

43 ■ 3947 

i(N) = 998093, g(10 6 ) = g(10 6 - 1) = - iV. 



n = 10 9 , = 2 14 3 9 5 6 7 5 ll 4 13 4 [17-31] 3 [37-263] 2 [269-150989] 

£(JV) = 999969437, 5 (10 9 ) = 5 (10 9 - 1) = 37 \ 150991 ^ 

v ; , yv ; ja ; 2 • 3 ■ 148399 



?i = 10 12 , N = 2 18 3 12 5 8 7 6 ll 5 13 5 [17-31] 4 [37-113] 3 [127-1613] 2 [1619-5476469] 

m = 999997526071, ,(10-) = 1621 ■ 1627 ■ 1637 ■ 5476483 

y ' ' y\ J 5475739-5476469 



n = 10 15 , N = 2 23 3 15 5 10 7 8 ll 7 13 6 17 6 [19-31] 5 [37-79] 4 [83-389] 3 

x[397-9623] 2 [9629- 192678817], 

£{N) = 999999940824564, 

, 192678823- 192678853- 192678883- 192678917 

g(10 i:5 ) = mlO - 1) = N. 

yv ; yv 7 389 • 9539 • 9587 • 9601 • 9619 • 9623 • 192665881 



Figure 6: The values g(n) for n = 10 6 , 10 9 , 10 12 , 10 15 . 



11 Open problems 

11.1 An effective bound for the benefit 

Let us define ben g(n) by (|6~T|) with AT and p defined by l|5TTjl and flUTO]). Is it 
possible to get an effective form of (|6.7j) . i.e., 

bengr(n) + n — £(g(n)) < Cp 

for some absolute constant C to determine? 

A hint is to apply Proposition [2] with M = f° r some r, where the 

Pi's are the r smallest primes not dividing N and the q^s are the 2r largest 
primes such that v qi (N) = 2, and, further, to apply effective results on the 
Prime Number Theorem like those of [28] or [5] . 

11.2 Increasing subsequences of g(n) 

An increasing subsequence of g is a set of k consecutive integers {n, n+l, . . . ,n+ 
k — 1} such that 

g(n - 1) = g(n) <g(n + l) < ... < g(n + k - 1) = g(n + fc). (11.1) 
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Due to a parity phenomenon!, these maximal sequences are rare. For n < 10 , 
there are only 9 values on n with k > 7. The record is n — 35464 with k = 20. 

Are there arbitrarily long maximal sequences? It seems to be a very difficult 
question. In [21] , (1.7), it is conjectured that there are infinitely many maximal 
sequences with k > 2. 

11.3 The second minimum 

Let us write gi(n) = g{n) > gs{n) > . . . > gi(n) = 1 all the integers such that, 
if a € & n , the order of a is equal to gi(n) for some i S {1,2,..., /}. From (|1.5jl . 
I is equal to the number of positive integers M satisfying i(M) < n. 

We might be interested in the computation of g-zin) or more generally, in 
the computation of g%{n) for 1 < i < iq where io is some (small) fixed constant. 

The basic algorithm (see Section [5} can be easily adapted for this purpose. 
It seems reasonnable to think that our algorithm, as sketched in II. 3[ can also 
be extended to get gi(n). 

11.4 Computing h(n) 

Let h(n) be the maximal product of primes pi x ,pi 2 , . . ■ ,Pi r under the condition 
Pii +Pi 2 + ■ • • + P» r <n(ris not fixed); h(n) can be interpreted as the maximal 
order of a permutation of the symmetric group & n such that the lengths of its 
cycles are all primes. 

A formula similar to (|1.2j) can be written: 

h(n) = max M. 

M squarefree 
l(M)<n 

The superchampion numbers are the product of the first primes. 

A related problem is to find an algorithm to compute h(n) for n up to 10 15 . 

11.5 Maximum order in GL(n, Z) 

Let G{n) be the maximum order of torsion elements in GL(n,Z). It has been 
shown in [TU] that 

Gin) = max M (11.2) 

L(M)<n 

where L is the additive function defined by L(l) = L(2) = and L(p a ) = 
ip(p a ) = p a -p"- 1 if p a > 3. 

From (TTl~2]l and (TO|) . it follows that g(n) < G(n) holds for all n's and it 
has been shown in |22j that limn—joo G(n) / g(n) = oo. 

Is it possible to adapt the algorithm described in this paper to compute G(n) 
up to 10 15 ? 
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